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Abstract:  We  propose  an  optimization-based  theoretical  approach  to 

tailor  the  optical  response  of  silver/silica  multilayer  nanospheres  over 
the  visible  spectrum.  We  show  that  the  structure  that  provides  the  largest 
cross-section  per  volume/mass,  averaged  over  a  wide  frequency  range,  is 
the  silver  coated  silica  sphere.  We  also  show  how  properly  chosen  mixture 
of  several  species  of  different  nanospheres  can  have  an  even  larger  minimal 
cross-section  per  volume/mass  over  the  entire  visible  spectrum. 
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1.  Introduction 

Nanoparticles  with  strong  optical  response,  characterized  by  scattering,  absorption  and  total 
cross-sections,  have  wide  applications  in  biomedical  imaging,  photothermal  therapy,  and  opti¬ 
cal  obscurance  [1-3].  Different  applications  require  different  optical  response  properties.  For 
instance,  real-time  biomedical  imaging  is  based  on  large  scattering  cross-sections,  while  pho¬ 
tothermal  therapy  requires  nanoparticles  with  large  absorption  cross-sections  and  small  scatte¬ 
ring  cross-sections.  For  obscurance  applications  [4, 5],  the  ideal  nanoparticle  should  typically 
have  large  total  cross-sections  over  the  whole  visible  spectrum  while  keeping  the  volume  or 
mass  of  nanoparticles  as  small  as  possible.  The  diversity  and  complexity  of  these  requirements 
necessitate  an  engineering  approach  of  nanoparticle  design. 

Previous  studies  on  the  optical  response  of  nanoparticles  are  mainly  based  on  parametric  ap¬ 
proach  [6-8],  which  works  well  for  simple  structures.  However,  when  the  structure  becomes 
complicated  and  the  number  of  design  parameters  increases,  optimization  becomes  the  prefer¬ 
able  approach  because  it  can  efficiently  explore  the  whole  parameter  space.  Furthermore,  be¬ 
cause  the  optimization  objective  function  can  be  an  arbitrary  transformation  of  the  frequency- 
dependent  cross-sections,  this  approach  is  very  powerful  in  tailoring  the  broadband  optical 
response  of  nanoparticles.  For  example,  some  applications  may  require  that  the  optical  reso¬ 
nance  has  both  a  strong  peak  value  and  a  wide  bandwidth.  This  can  be  achieved  by  maximizing 
the  average  cross-section  over  the  bandwidth  of  interest.  For  obscurance  applications,  we  want 
the  total  cross-section  to  be  consistently  large  over  the  whole  visible  spectrum.  This  is  equiva¬ 
lent  to  maximizing  the  minimal  cross-section  over  this  spectrum.  In  this  paper,  we  will  use  an 
optimization  tool  to  tailor  the  optical  response  of  multilayer  nanospheres  over  wide  frequency 
range  of  interest. 

Before  we  start,  we  need  to  select  the  material  system.  Nanoparticles  composed  of  metal 
and  dielectric  materials  support  surface  plasmons  on  the  metal/dielectric  interfaces  and  can 
strongly  interact  with  light  in  the  visible  range  [9-13].  At  resonance,  the  cross-sections  of 
these  nanoparticles  are  much  larger  than  their  physical  cross-sections,  which  makes  them  super- 
scatterers  and  super-absorbers  [14-17].  Furthermore,  the  plasmon  resonance  frequency  can  be 
tuned  by  varying  the  physical  structure  of  the  nanoparticles.  In  order  to  be  able  to  tailor  the 
optical  response  of  nanoparticles  over  a  wide  frequency  range,  we  choose  the  metal/dielectric 
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material  system.  For  concreteness,  we  will  focus  on  the  silver/silica  material  system. 

In  this  article,  we  first  formulate  the  Transfer  Matrix  Method  to  calculate  the  optical  response 
of  multilayer  nanospheres.  After  briefly  reviewing  the  optical  response  of  bilayer  silver/silica 
nanospheres,  we  proceed  to  optimize  the  average  cross-section  of  various  silver/silica  multi¬ 
layer  nanospheres.  Our  optimization  results  show  that  the  structure  with  the  maximal  average 
cross-section  is  the  bilayer  silver/silica  structure  with  silver  as  the  shell.  Finally,  we  investigate 
using  a  mixture  of  several  species  of  bilayer  nanospheres  to  enhance  the  minimal  cross-sections 
over  the  entire  visible  range. 

2.  Calculation  of  optical  response  via  transfer  matrix  method 


Fig.  1.  Schematic  of  an  n  layer  nanosphere  embedded  in  infinite  dielectric  medium.  The 
outer  radius  and  dielectric  function  of  individual  layers  are  (/?/,£;),/  =  1,2 ,  The  di¬ 

electric  function  of  the  medium  is  £m.  The  solid  lines  represent  an  incident  plane  wave 
which  contains  incoming  and  outgoing  waves.  The  dashed  line  represents  the  scattered 
wave  which  only  contains  outgoing  wave. 

We  start  by  formulating  the  Transfer  Matrix  Method  for  multilayer  nanosphere  [18, 19].  Con¬ 
sider  a  multilayer  nanosphere  shown  in  Fig.  1 .  Because  of  the  spherical  symmetry,  the  fields  at  a 
given  incident  frequency  can  be  decomposed  into  two  orthogonal  polarizations:  transverse  elec¬ 
tric  (TE)  and  transverse  magnetic  (TM).  For  TE  polarization,  the  electric  fields  can  be  written  as 
E te  =  V  x  y(\>te  •  For  TM  polarization,  the  magnetic  fields  can  be  written  as  H tm  =  V  x  y(\>tm- 
The  scalar  potential  (\>te  and  satisfy  the  scalar  Helmholtz  equation  V20  +  &20  =0  where 
k 2  =  £02£(r).  Due  to  the  spherical  symmetry,  0  can  be  decomposed  into  a  discrete  set  of  spher¬ 
ical  modes:  0/m  =  Ri(r)pjm^  (cos  6)  exp (imcp)  with  /  =  0,1,2, ...  and  m  =  Since  £(r) 

is  a  constant  £*•  inside  the  ith  shell,  Ri(r)  is  a  linear  combination  of  the  first  and  second  kind 
spherical  Bessel  functions  within  the  individual  shells: 


Ri(r)\i=AiJi(kir)  +5iy/(^r) 


(1) 


The  coefficients  (A^Bi)  of  adjacent  shells  are  linked  by  the  transfer  matrix  of  the  interface: 


^i+ 1 
Bi+ 1 


(2) 


The  matrix  element  is  determined  by  the  boundary  condition  satisfied  by  Ri(r),  which  comes 
from  the  continuity  of  the  tangent  components  of  E  and  H  across  the  boundary.  For  TE  po- 
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larization,  r/?;(r)  and  ( rRi{r ))'  are  continuous  across  the  boundary.  By  writing  the  continuity 
conditions  in  matrix  form,  we  get: 


Mi+ i,j 


jl(ki+lRi)  yi(ki+iRi) 

j'l(ki+\Ri)ki+\Ri  +  ji(ki+\Ri)  ^(ki+tfiiki+iRi+yiiki+iRi) 

jlikiRi)  yi{kiRi) 

jl  ( kiRi)kiRi  +  ji  (hRi)  y\ {klRi)klRl  +  yt  (*,■/?,■) 


For  TM  polarization,  rRi(r)  and  (e  1  rRi(r))'  are  continuous  across  the  boundary.  By  writing 
the  continuity  conditions  in  matrix  form,  we  get: 


Mi+hi 


jiiki+iRi) 

/t(ki+lRi)ki+lRi  +  ji(ki+lRi) 


yi(ki+iRi) 

y'i(ki+iRi)ki+iRi+yi(ki+lRi) 


jiikiRt) 

^UifcRdkiRi  +  MkiRi)) 


yiikiRi) 

^WWdWi+yifaRi)) 


(4) 


The  transfer  matrix  of  the  whole  system  can  be  calculated  by  cascading  the  transfer  matrices  of 
individual  interfaces. 


An+i 
Br+ l 


—  Mn-\-\,nMnn—i  •••^3,2^2,! 


Ai 

Bi 


(5) 


Since  the  second  kind  of  Bessel  function  is  singular  at  the  origin,  we  can  set  A\  =  1  and 
B\  =0.  So  the  coefficients  of  Bessel  functions  in  the  surrounding  medium  are  directly  given 
by  the  transfer  matrix  element,  An+\  =  M\\  and  Bn+ 1  =  M21.  Within  the  surrounding  medium, 
it  is  convenient  to  write  the  radical  function  as  a  linear  combination  of  the  spherical  Hankel 
functions: 

Rl(r)  |„+i  =  Cn+ihj(kn+lr)  +Dn+1hf(kn+1r)  (6) 

Taking  the  convention  that  the  fields  vary  in  time  as  e~l(0t ,  hj  (kn+\r)  and  hj(kn+\r)  correspond 
to  outgoing  and  incoming  waves  respectively.  The  reflection  coefficient  of  the  whole  system  is 
given  by 

Cn+ 1  M11—M21 

ri  = - = -  (7) 

Dn+1  Mn  +  iM2\ 

The  reflection  coefficient  as  a  function  of  frequency  determines  the  optical  response  of  the 
nanoparticle  under  all  possible  illumination  conditions.  Specifically,  when  the  nanoparticle  is 
illuminated  by  a  linearly  polarized  plane  wave,  the  incident  field  can  be  decomposed  into  both 
TE  and  TM  channels  (/,ra)  with  /  =  1,2,...  and  m  —  — 1,1.  For  each  channel,  the  incident  field 
contains  both  incoming  and  outgoing  waves  carrying  the  same  power  [14-17] 

Pt,m=±\  =  ^(2/ +l)/0  (8) 

where  To  is  the  incident  intensity,  and  A  is  the  wavelength  in  the  surrounding  medium.  The 
scattered  field  is  a  purely  outgoing  wave  characterized  by  the  scattering  coefficient  S\  which  is 
related  to  the  reflection  coefficient  through  Si  =  (77  —  l)/2.  The  scattered  and  absorbed  power 
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in  this  channel  is  given  by: 


P?:=±  i  =  ^(2/+l)/0|l-r;|2  (9) 

PftL±  i  =  ^(2/+l)/0(l-|n|2)  (10) 

By  summing  over  the  contribution  from  all  channels  of  TE  and  TM  polarization,  we  get  the 
scattering  and  absorption  cross-sections: 

c>sca  =  XXir(2/+1)l1-?wl2  (11) 

a  /=1  871 

°abs  =  l£^(2/+l)(l-k<T,/|2)  (12) 

where  cj  is  TE  or  TM.  The  total  cross-section  is  the  sum  of  the  scattering  and  absorption  cross- 
sections.  CJ \ot  —  GSca  H-  Gabsm 

3.  Optical  response  of  silver/silica  bilayer  nanospheres 


A  (nm) 


<3 

a 

<X> 

Oh 


o 

■+o 

b 


Fig.  2.  The  total  cross-section  of  silica  coated  silver  spheres  suspended  in  ethanol.  The 
cross-section  is  normalized  by  volume  (the  left  axis)  and  mass  (the  right  axis).  The  insert 
is  a  TEM  image  of  the  fabricated  nanoparticles.  The  radius  of  the  silver  core  has  a  distribu¬ 
tion  with  mean  26.3nm  and  standard  deviation  9.3nm.  The  thickness  of  the  silica  shell  is 
around  25.3nm.  The  red  line  is  the  measured  total  cross-section.  The  black  bar  represents 
the  standard  deviations  from  eight  transmission  measurement  on  eight  samples.  The  blue 
line  is  the  Transfer  Matrix  calculation  of  the  total  cross-section  with  the  radius  of  the  silver 
core  sampled  from  the  measured  distribution  and  the  thickness  of  the  silica  shell  fixed  at 
25.3nm.  The  dielectric  function  of  ethanol  is  taken  as  em  =  1.85. 
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In  this  section,  we  analyze  the  optical  response  of  silver/silica  bilayer  nanospheres  as  the  build¬ 
ing  elements  of  multilayer  structures.  For  silver,  the  complex  dielectric  function  as  a  function 
of  frequency  is  generated  by  linearly  interpolating  the  experimental  data  [20] .  For  simplicity, 
the  size  dependence  of  Ag’s  dielectric  function  is  not  taken  into  account.  For  silica,  the  dielec¬ 
tric  function  is  taken  as  a  constant  £  =  2.1.  There  are  two  configurations  of  silver/silica  bilayer 
nanoparticle  depending  on  the  core  material. 

First,  we  consider  silica  coated  silver  spheres.  Colloidal  suspensions  of  such  nanoparticles 
were  obtained  from  Nanocomposix  (Nanocomposix  Inc.,  San  Diego,  CA).  Figure  2  shows  the 
measured  and  calculated  total  cross-sections  of  the  fabricated  nanoparticles.  The  calculation 
agrees  quite  well  with  the  measurement.  The  total  cross-section  peaks  around  455nm.  This  peak 
comes  from  the  /  =  1  surface  plasmon  mode  at  the  silver/silica  boundary.  The  peak  wavelength 
only  varies  slightly  when  the  inner  and  outer  radius  change.  For  instance,  consider  a  silica 
coated  silver  sphere  suspended  in  air.  Fixing  the  outer  radius  at  50nm,  the  peak  wavelength 
varies  from  410nm  to  415nm  when  the  inner  radius  varies  from  5nm  to  45nm.  Fixing  the  aspect 
ratio  R1/R2  at  0.8,  the  peak  wavelength  varies  from  390nm  to  480nm  when  the  outer  radius 
varies  from  25nm  to  75nm. 

In  contrast,  the  surface  plasmon  resonance  of  the  reverse  configuration  has  great  tunability 
over  the  visible  range  [9-13].  Consider  a  silver  coated  silica  sphere  suspended  in  air.  Fixing  the 
outer  radius  at  50nm,  the  peak  wavelength  varies  from  405nm  to  720nm  when  the  inner  radius 
varies  from  5nm  to  45nm.  Besides  the  peak  wavelength,  the  relative  strength  of  scattering  and 
absorption  cross-sections  in  the  total  cross-section  also  vary.  For  [R\,R2\  =  [5nm,  50nm],  the 
absorption  cross-section  accounts  for  25%  of  the  total  cross-section  at  resonance.  For  [R\ ^2]  = 
[45nm,  50nm],  this  percentage  rises  to  60%.  The  tunability  of  the  resonance  wavelength  and 
the  tunability  of  the  total  cross-section  composition  makes  silver  coated  silica  sphere  a  good 
candidate  for  achieving  broadband  optical  response. 

4.  Optimization  of  average  cross-sections  over  wide  frequency  range 

The  optical  response  of  silver/silica  bilayer  nanosphere  indicates  that  increasing  the  number  of 
metal/dielectric  interfaces  can  provide  additional  tunability  in  the  optical  response.  Using  this 
insight,  we  aim  to  design  silver/silica  multilayer  nanosphere  with  a  large  average  cross-section 
over  wide  frequency  range.  The  figure  of  merit  (FOM)  is  the  scattering,  absorption,  and  total 
cross-section  averaged  over  the  target  frequency  range,  normalized  by  volume  or  mass. 

]  r  (timax 

FOM  =  —  /  GnormalizeddCfi  (12) 

(Umax  —  uJmin  J  (omin 

For  concreteness,  we  take  400-600nm  and  600-800nm  as  the  target  frequency  range  of  interest. 
The  structure  under  consideration  is  a  multilayer  nanosphere  with  alternating  silver  and  silica 
layers  (six  layers  in  total).  The  design  parameters  are  the  thicknesses  of  individual  layers.  The 
lower  bounds  of  the  thicknesses  are  set  to  be  zero.  The  upper  bound  of  the  allowed  thickness 
is  set  to  be  a  large  value  (1/im).  Therefore,  this  general  structure  includes  structures  with  fewer 
layers  (one  through  five  layers)  as  boundary  points. 

We  performed  the  optimization  using  numerical  optimization  package  NLOPT  [21].  Since 
this  problem  is  nonconvex,  there  are  many  local  optima.  To  find  the  global  optima  in  the  design 
parameter  space,  we  used  Multi-Level- Single-Linkage  (MLSL)  algorithm.  This  algorithm  per¬ 
forms  a  sequence  of  local  optimization  from  random  points  by  a  clustering  heuristic  that  helps 
it  to  avoid  repeated  searches  of  the  same  optima  [22] .  The  local  optimization  algorithm  used 
here  is  BOBYQA  [23].  This  algorithm  performs  derivative-free  bound-constrained  optimiza¬ 
tion  using  an  iteratively  constructed  quadratic  approximation  of  the  objective  function. 

Figure  3  summarizes  the  optimization  results.  In  all  cases,  the  optimal  structure  returned  by 
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Fig.  3.  Optimization  of  average  cross-sections  over  wide  frequency  range.  The  structure 
under  consideration  is  a  silver/silica  multilayer  nanosphere.  The  optimal  structure  found 
by  the  optimization  engine  is  always  silver  coated  silica  sphere.  For  all  subfigures,  blue 
(red)  lines  show  the  optimized  average  cross-sections  over  the  blue  (red)  shaded  frequency 
range.  (a)(b)(c)  correspond  to  scattering,  absorption  and  total  cross-sections  per  volume 
respectively.  (d)(e)(f)  correspond  to  scattering,  absorption  and  total  cross-sections  per  mass 
respectively.  The  radius  of  the  silica  cores  and  the  thickness  of  silver  shells  exhibiting  the 
cross-sections  shown  above  are  given  in  Table  1. 


the  optimization  engine  is  always  a  silver  coated  silica  sphere.  Although  multilayer  structures 
can  offer  greater  tunability  of  the  optical  response,  bilayer  structures  already  maximize  the 
average  cross-section  over  wide  frequency  range.  From  Table  1,  we  can  see  that  nanospheres 
with  an  outer  radius  around  70nm  have  the  largest  normalized  average  scattering  cross-section. 
The  wavelength  of  the  scattering  peak  can  be  further  tuned  by  varying  the  aspect  ratio.  For  ab¬ 
sorption  cross-sections,  our  optimization  engine  found  many  local  optima  with  approximately 
equal  FOM’s.  These  local  optima  have  the  same  aspect  ratio,  and  the  thickness  of  the  silver 
layer  varies  from  zero  to  several  nanometers.  This  can  be  explained  by  the  quasi-static  approxi¬ 
mation.  The  absorption  cross-section  of  a  nanoparticle  can  be  written  as  oabs(co)  =  coIm[a(co )], 
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where  a (co)  is  the  polarizability  [24-26].  When  the  nanoparticle  diameter  is  much  smaller  than 
the  wavelength,  the  quasi-static  approximation  is  valid,  under  which  the  polarizability  is  pro¬ 
portional  to  the  volume  of  the  nanoparticle  with  the  proportionality  coefficient  dependent  on 
frequency  and  aspect  ratio  [24].  Therefore,  the  normalized  absorption  cross-section  averaged 
over  a  frequency  range  is  only  determined  by  the  aspect  ratio  and  independent  of  the  nanopar¬ 
ticle  diameter  as  long  as  the  nanoparticle  diameter  is  much  smaller  than  the  wavelength.  In 
Table  2,  the  thickness  of  the  silver  layer  is  set  to  be  >  2nm.  Nanospheres  with  thinner  silver 
layers  have  approximately  the  same  normalized  average  absorption  cross-section  but  are  dif¬ 
ficult  to  fabricate.  The  structures  that  give  the  largest  average  scattering  and  absorption  cross- 
sections  are  quite  different.  Since  the  FOM  of  the  absorption  cross-section  is  about  twice  as 
large  as  the  FOM  of  the  scattering  cross-section,  the  structure  that  gives  the  largest  average 
total  cross-section  is  essentially  identical  to  the  structure  providing  the  largest  average  absorp¬ 
tion  cross-section.  The  structures  given  in  Table  1  are  super-scatters  and  super- absorbers.  For 
instance,  the  structure  with  [R,T]  =  [60.40nm,  8.68nm]  has  an  average  scattering  cross-section 
of  8.65m2/g  over  600-800nm,  while  its  physical  cross-section  is  only  2.07m2/g.  The  structure 
with  [/?,  T]  =  [18.09nm,  2.00nm]  provides  an  average  total  cross-section  of  17.52m2/g  over 
600-800nm,  which  means  that  only  lg  of  such  nanoparticles,  when  fully  dispersed,  can  ob¬ 
scure  an  area  as  large  as  17.52m2. 


Table  1.  Optimization  of  average  cross-sections 


Normalized  by  volume 


cross-section 

range  (nm) 

silica  (nm) 

silver  (nm) 

FOM  (1/nm) 

Gsca 

400-600 

31.25 

26.65 

0.0486 

600-800 

60.32 

9.65 

0.0464 

Gabs 

400-600 

6.07 

2.00 

0.0767 

600-800 

14.80 

2.00 

0.0817 

Gfot 

400-600 

6.09 

2.00 

0.0773 

600-800 

16.12 

2.00 

0.0846 

Normalized  by  mass 


cross-section 

range  (nm) 

silica  (nm) 

silver  (nm) 

FOM  (m2/g) 

Gsca 

400-600 

38.53 

17.88 

5.64 

600-800 

60.40 

8.68 

8.65 

Gabs 

400-600 

6.54 

2.00 

10.87 

600-800 

17.89 

2.00 

16.71 

Gtot 

400-600 

8.71 

2.00 

10.93 

600-800 

18.09 

2.00 

17.52 

5.  Optimization  of  the  minimal  cross-sections  over  wide  frequency  range 

In  this  section,  we  aim  to  design  nanoparticles  with  consistently  large  cross-sections  over  a  wide 
frequency  range.  This  is  equivalent  to  maximizing  the  minimal  cross-section  over  the  target 
frequency  range.  From  the  previous  section,  we  see  that  silver  coated  silica  spheres  have  strong 
surface  plasmon  resonances,  and  the  peaks  of  their  resonances  are  highly  tunable.  Therefore,  we 
expect  that  a  mixture  of  several  species  of  silver  coated  silica  spheres  with  different  resonance 
frequencies  can  effectively  cover  a  wide  frequency  range. 

This  intuition  can  be  formalized  by  optimization  language.  We  take  the  target  wavelength 
range  to  be  400-800nm.  The  FOM  is  the  minimal  cross-section  over  this  range,  normalized  by 
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Fig.  4.  Optimization  of  minimal  cross-sections  over  a  wide  frequency  range.  The  structure 
under  consideration  is  a  mixture  of  several  species  of  silver  coated  silica  spheres.  The  target 
frequency  range  is  shaded  in  yellow.  For  all  subfigures,  blue,  red,  black  lines  corresponds  to 
one,  two,  and  three  species  of  nanospheres.  The  black  dashed  lines  in  (c)  and  (f)  correspond 
to  ten  species  of  nanospheres.  (a)(b)(c)  correspond  to  scattering,  absorption,  and  total  cross- 
sections  per  volume  respectively.  (d)(e)(f)  correspond  to  scattering,  absorption  and  total 
cross-sections  per  mass  respectively.  The  radius  of  the  silica  cores  and  the  thickness  of 
silver  shells  corresponding  to  these  cross-sections  are  given  in  Table  2. 

either  volume  or  mass. 

FOM  =  min  Onormanzeci  (14) 

CO 

The  structure  under  consideration  is  a  mixture  of  N  species  of  silver  coated  silica  spheres  where 
N  =  1,2,3,....  The  design  parameters  are  the  size  parameters  of  individual  species  and  relative 
weights  (i.e.  proportions)  of  each  species  in  the  mixture.  The  weights  represent  the  relative 
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weights  in  volume  (mass)  when  the  normalization  is  over  volume  (mass). 

N 

O normalized  x  ,  W jOj, normalized  (1^) 

i=  1 

When  the  size  parameters  of  individual  species  are  fixed,  the  problem  of  finding  the  optimal 
weights  turns  out  to  be  a  Linear  Programming  (LP)  program.  Therefore,  we  employed  a  two- 
level  optimization  structure.  In  the  lower  level,  we  used  a  standard  LP  solver  to  find  out  the 
optimal  weights  given  the  current  size  parameters.  The  resulting  FOM  as  a  function  of  size 
parameters  is  further  optimized  in  the  upper  level  with  the  same  nonlinear  method  we  used  in 
the  previous  section.  This  separation  into  linear  and  nonlinear  parts  of  the  original  optimization 
problem  reduces  the  dimension  of  the  parameter  space  and  helps  MLSL  algorithm  to  find  the 
global  optimal  in  less  iterations. 

Figure  4  summarizes  the  results.  We  can  see  that  the  optimization  engine  tries  to  build 
a  plateau  over  the  target  range  to  maximize  the  minimal  cross-section.  For  scattering  cross- 
sections,  it  is  relatively  easy  to  build  such  a  plateau  since  the  scattering  peaks  have  wide  band¬ 
width.  When  there  is  only  one  species  (the  blue  lines  in  Fig.  4(a)(d)),  the  nanoparticle  has  a 
large  size  so  that  its  /  =  1  resonance  can  cover  the  large  wavelength  region  and  its  7  =  2  reso¬ 
nance  can  cover  the  small  wavelength  region.  When  the  number  of  species  increases,  the  new 
species  try  to  cover  the  dips  in  the  original  scattering  spectra  with  their  resonant  peaks.  On 
the  other  hand,  it  is  relatively  difficult  to  build  an  absorption  plateau  because  absorption  peaks 
have  narrow  bandwidth.  The  FOM  of  the  scattering  cross-section  is  about  twice  as  large  as  that 
of  the  absorption  cross-section.  The  total  cross-section  can  be  enhanced  significantly  (35%  for 
volume  normalization  and  46%  for  mass  normalization)  when  N  increases  from  1  to  2.  The  en¬ 
hancement  when  N  increases  from  2  to  3  is  only  moderate.  The  benefit  of  adding  more  species 
gradually  saturates. 

6.  Concluding  remarks 

In  this  article,  we  used  optimization  tools  to  tailor  the  optical  response  of  silver/silica  multilayer 
spheres.  We  show  that  the  structure  that  gives  the  largest  average  cross-section  over  wide  fre¬ 
quency  range  is  the  bilayer  structure  with  a  silver  shell.  We  also  show  that  using  several  species 
of  nanoparticles  can  significantly  enhance  the  minimal  cross-section  over  the  whole  visible 
range  although  this  enhancement  saturates  when  the  number  of  species  increase.  Because  the 
FOM  can  be  an  arbitrary  function  of  the  frequency  dependent  cross-sections,  the  optimization 
approach  described  here  can  be  used  to  design  nanoparticles  with  more  complicated  optical 
response.  We  are  investigating  efficient  numerical  methods  of  calculating  and  optimizing  the 
optical  response  of  non- spherical  geometries  as  future  directions. 

This  research  was  funded  by  the  Department  of  the  Army  Basic  Research  Program  and  spon¬ 
sored  by  the  Edgewood  Chemical  Biological  Center.  Support  was  also  provided  by  the  U.S. 
Army  Research  Office  under  contract  W911NF-07-D-004. 
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Table  2.  Optimization  of  minimal  cross-sections 


Normalized  by  volume 


cross-section 

#  of  species 

silica  (nm) 

silver  (nm) 

weights  (%) 

FOM  (1/nm) 

Gsca 

1 

68.09 

55.44 

100 

0.0192 

2 

52.54 

44.21 

82.8 

0.0229 

73.41 

7.63 

17.2 

3 

52.85 

47.27 

74.4 

0.0235 

75.61 

8.31 

20.8 

35.55 

23.21 

4.8 

Gabs 

1 

60.14 

10.00 

100 

0.0032 

2 

37.03 

4.00 

57.5 

0.0075 

44.02 

15.99 

42.5 

3 

22.95 

10.26 

37.4 

0.0114 

26.08 

4.39 

34.3 

21.96 

2.00 

28.3 

Gtot 

1 

80.15 

41.87 

100 

0.0201 

2 

59.96 

34.12 

82.8 

0.0270 

61.96 

6.29 

17.2 

3 

50.21 

32.98 

75.2 

0.0285 

52.97 

4.90 

14.4 

55.76 

7.32 

10.4 

10 

0.0310 

Normalized  by  mass 


cross-section 

#  of  species 

silica  (nm) 

silver  (nm) 

weights  (%) 

FOM  (m2/g) 

Gsca 

1 

67.95 

52.38 

100 

2.10 

2 

52.77 

38.70 

87.4 

2.80 

75.72 

8.17 

12.6 

3 

76.71 

21.93 

46.1 

2.88 

21.96 

34.54 

43.2 

70.45 

7.31 

10.7 

Gabs 

1 

42.80 

6.55 

100 

0.58 

2 

20.56 

2.15 

51.3 

1.31 

37.16 

13.61 

48.7 

3 

27.43 

14.22 

36.9 

1.59 

49.06 

11.02 

34.9 

36.66 

3.66 

28.2 

Gtot 

1 

92.95 

39.04 

100 

2.46 

2 

55.54 

29.60 

83.5 

3.59 

67.77 

7.39 

16.5 

3 

53.36 

29.01 

60.1 

3.79 

79.49 

10.25 

31.1 

24.64 

25.04 

8.8 

10 

4.62 
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